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Abstract. We present the analysis of a new flux-limited 
sample of bright quasars and Seyfert 1 galaxies in an ef- 
fective area of 611 deg 2 , drawn from the Hamburg/ESO 
survey. We confirm recent claims that bright quasars have 
a higher surface density than previously thought. Special 
care was taken to avoid morphological and photometric 
biases against low-redshift quasars, and about 50% of 
the sample objects are at z < 0.3, spanning a range of 
AM ~ 8 in absolute magnitudes. While our derived space 
densities for low-luminosity Seyfert 1 nuclei are consistent 
with those found in the literature, we find that luminous 
QSOs, Mb < —24, are much more numerous in the lo- 
cal universe than previous surveys indicated. The opti- 
cal luminosity functions of Seyfert 1 nuclei and QSOs join 
smoothly, and if the host galaxy contributions are taken 
into account, a single power-law of slope a = —2.2 de- 
scribes the combined local luminosity function adequately, 
over the full range in absolute magnitude. Comparing our 
data with published results at higher redshifts, we can rule 
out pure luminosity evolution as an acceptable parametri- 
sation; the luminosity function of quasars changes shape 
and slope with z, in the sense that the most luminous 
quasars show the weakest evolution. 

Key words: Galaxies: luminosity function - Quasars: 
general - Galaxies: Seyfert 



1. Introduction 

The average luminosities of quasars, and possibly many 
other properties, depend strongly on cosmological epoch. 
Shape, normalisation, and evolution of the quasar lumi- 
nosity function (QLF) are among the most basic descrip- 
tions of the quasar population. A good knowledge of the 
local (z w 0) QLF is essential for at least two purposes: 
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first, to serve as zero-point for quasar evolution studies; 
second, to provide a reference distribution law for unbi- 
ased subsamples constructed for statistical investigations 
(e.g., for spectral properties, host galaxies, etc). 

Unfortunately, current determinations of the QLF are 
limited mainly to z ^ 0.3, for the simple reason of lack- 
ing appropriate survey data. Most optical QSO surveys 
discriminate against objects with extended morphology, 
causing severe incompleteness (and bias) at low redshifts. 
At the low-luminosity end of the QSO-Seyfert family, 
dedicated galaxy surveys have prompted investigations of 
the Seyfert 1 luminosity function (e.g., Cheng et al. 1985; 
Huchra & Burg 1992). However, these surveys are neces- 
sarily incomplete at high nuclear luminosities, where the 
active nuclei outshine the surrounding hosts, and the ob- 
jects are not selected as galaxies anymore. 

In 1990 we started the 'Hamburg/ESO survey' (HES), 
a new wide-angle survey for bright QSOs and Seyferts, de- 
signed to avoid the above selection biases as far as possible. 
A detailed description of the HES was given by Wisotzki 
et al. (1996; hereafter Paper I), and a first list of 160 
newly discovered QSOs was published by Reimers et al. 
(1996; Paper II). Basically, QSO candidates are selected 
by automated procedures from digitised objective-prism 
plates taken with the ESO Schmidt telescope, with subse- 
quent slit spectroscopy of all candidates above the mag- 
nitude limit of typically B\ lm ~ 17.0-17.5, over an area 
of ~ 5000 deg 2 at high galactic latitude. To prevent mor- 
phological bias against sources with resolved structure, 
such as Seyferts or gravitational lenses, objects with non- 
stellar appearance are not excluded from the candidate 
lists. Among other goals, we aim to construct a large, op- 
tically flux- limited sample of low-redshift 'Type 1 AGN', 
including QSO/Seyfert borderline objects. 

In this paper we publish the results of an analysis of 
a first part of the HES, based on 33 Schmidt fields, with 
an effective area of 611 deg 2 . We present the first direct 
determination of the local quasar luminosity function from 
a single optical survey, and discuss some implications of 
our results. A comprehensive description of the analysis 
has been given by Kohler (?). We adopt a cosmological 
model with Hq = SOkms" 1 Mpc^ 1 , qo = 0.5, and A = 0. 
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Table 1. Properties of surveyed fields. Field designations are 
given as standard ESO/SERC field numbers. Nh is the col- 
umn density of Galactic neutral hydrogen in the field centre, in 
fO 20 cm~ 2 . Ab is the corresponding extinction in the B band. 
Q e ff is the effective survey area for the field, in deg 2 . 
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2. The quasar sample 

2.1. Survey area 



The investigated area was denned by 33 ESO Schmidt 
fields located in the region -28° < S < -7° and 9 h 40 m < 
a < 14 h 40 m , thus all in the North Galactic hemisphere. 
A list of the survey fields and their properties is given in 
Table 1. 

In each field, one direct plate from the ESO (B) atlas, 
and one HES spectral (objective prism) plate were avail- 
able. The formal total area subtended by 33 ESO plates is 
~ 25 deg 2 x 33 = 825 deg 2 . Losses of usable area occurred 
because (1) plates overlapped in adjacent fields, especially 
if not precisely centred; (2) occasionally direct and spec- 
tral plates were positionally mismatched, and only the 
common area could be used; (3) most important, overlap- 
ping spectra (and to a lesser degree, direct images) ren- 
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Fig. 1. Effective area of the surveyed region as a function of 
B magnitude (without correction for Galactic extinction). 



dered a certain percentage of spectra unprocessable. These 
losses have been quantified and incorporated into 'effec- 
tive areas' fl e g for each field, given in Table 1. Because 
of the field-to-field variations in limiting magnitude, the 
total effective survey area is a function of apparent magni- 
tude, as plotted in Fig. 1. The maximum area is 611 deg 2 
for B < 16.64 and decreases gradually to zero for sources 
fainter than B = 17.66. 

While many quasar surveys were carried out at very 
high Galactic latitudes where foreground extinction cor- 
rections are small (although systematically non-zero), this 
is not the case for our present set of fields. We have used 
the data from Stark et al. (1992) to obtain the neutral 
hydrogen column density in the centre of each field, and 
converted this into a B band extinction using the formula 
A B = 4.2 xNh/59 where N H is given in units of 10 20 cm~ 2 
(cf. Spitzer 1978); Table 1 lists N H and A B for each field. 
The average extinction is 0.36 mag. 

2.2. Sample selection 

The limiting brightness of QSO candidates during the au- 
tomated selection procedure was initially defined by a 
minimum S/N ratio of ~ 3 in the digital spectra. We 
later found, as expected (cf. Paper I), that close to the 
selection limit the object lists became incomplete. Sev- 
eral experiments showed that systematic incompleteness 
occurred only for the faintest spectra, and that the effect 
disappeared for a limiting S/N > 5. We translated this 
S/N cutoff into a magnitude limit on the direct plate, us- 
ing the calibration of photographic magnitudes described 
below. For the 33 fields, the values for B\i m range between 
17.66 and 16.64, depending largely on the quality of the 
spectral plates. 

Altogether, 115 QSOs and Seyfert 1 galaxies were de- 
tected by the HES selection criteria in these fields (pri- 
marily different measures of UV excess; cf. Paper I) that 
were either confirmed by our own follow-up spectroscopy 
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Table 2. The flux-limited sample of QSOs with z > 0.07. 
The entries in column B marked by colons are photographic, 
all other magnitudes are CCD measurements corrected to the 
zero-points of the photographic plates. Column B lists the 
extinction-corrected magnitudes. 

Name z B Field B M B 



(Paper II) , or previously catalogued as QSO or Seyfert 1 
in the compilation of Veron-Cetty & Veron (1993). Note 
that in this paper we deliberately ignore the traditional 



subdivision Seyferts vs. 'real' QSOs based on some arbi- 
trary luminosity threshold. A lower redshift limit of 0.07 
was applied to avoid possible incompleteness due to host 
galaxy contamination (but see the z < 0.07 Seyfert sam- 
ple below). Several of the lowest redshift objects would 
nevertheless be called Seyferts by many others. 

There remained 48 QSOs with z > 0.07 located within 
the effective survey area and above the completeness lim- 
iting magnitude. These form the flux-limited sample listed 
in Table 2. A comparison with the Veron-Cetty & Veron 
(1993) catalogue revealed no additional QSOs that should 
have been included in the flux-limited sample. This test for 
'completeness' is, however, not very strong, as few other 
surveys have covered this region of the sky, and ~ 90 % of 
the sample objects are new discoveries made by the HES. 

2.3. Photometry 

Obtaining consistent and unbiased photometry of low- 
rcdshift quasars with detectable host galaxies is not triv- 
ial. During the survey phase, normally only photographic 
data, in our case from the digitised direct B plates, are 
available. These plates usually reach much fainter magni- 
tudes than objective-prism plates, but they also saturate 
earlier, especially if not the original plates but glass copies 
are used (as was the case for several of our fields). The 
conventionally employed isophotal or large aperture mag- 
nitudes, although well suited for point sources, tend to give 
much larger total luminosities for QSOs with detectable 
hosts than desired for an investigation of the nuclear lumi- 
nosity function. Additional complications can result from 
non-linear filtering effects inherent in the digitisation pro- 
cedure. 

For the present sample we have developed a two-step 
procedure to obtain unbiased and survey-consistent mag- 
nitudes. The photographic plates were calibrated by CCD 
sequences obtained with the ESO 90 cm telescope in 1993 
and 1994. These sequences will be published in a separate 
paper. The sample selection was based on simulated di- 
aphragm photometry with an aperture of the size of the 
seeing disk. For point sources, an accurate calibration re- 
lation with very little intrinsic scatter, ~ 0.15 mag, could 
be derived from the CCD sequences (which contained only 
stars). These small-aperture magnitudes correlated also 
well with the corresponding S /N ratios of the digital spec- 
tra, important for an accurate definition of limiting mag- 
nitudes. 

In a second step, we obtained CCD photometry for 
most of the sample QSOs. Exposures in the B band were 
available from the mandatory acquisition images at the 
ESO 3.6 m and 2.2 m telescopes (cf. Paper II). After stan- 
dard reduction and colour term corrections, the photomet- 
ric zero points of these images were shifted to the system of 
the corresponding photographic plate by using surround- 
ing unsaturated field stars as references. The resulting 
corrected CCD magnitudes, again measured in an aper- 
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Table 3. Seyfert 1 galaxies with z < 0.07. The CCD magnitudes are for the small aperture (see text). Absolute magnitudes are 
given as observed (M fl ^ s ), and with host galaxy contribution M ga[ subtracted (M nuc ). 
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ture of diameter approximately equal to that of the seeing 
disk, should be much more accurate than the photographic 
data, at least for extended objects, while magnitudes of 
point sources should be more or less unchanged. We tested 
this assertion by comparing photographic and zero-point 
corrected CCD photometry of the z > 0.3 QSOs (which 
should be true point sources) and found a mean differ- 
ence of only 0.01 mag. For the further analysis we used 
the zero-point corrected B band CCD magnitudes. For 
three objects no CCD data were available, and for those 
we substituted the photographic measurements (cf. Ta- 
ble 2). 



2.4- The additional Seyfert sample 

To obtain further insight into the continuity of properties 
between high-luminosity QSOs and low-luminosity Seyfert 
nuclei, we conducted a dedicated subsurvey for Seyfert 
galaxies, targeting both type 1 and type 2 objects. The 
criteria were different from those used in the main sur- 
vey, adapted to discriminate between active and inactive 
galaxies rather than between stars and QSOs. More de- 
tails about this subsurvey will be given in a later paper 
(Kohler et al., in prep.) where we shall also present our 
results concerning Seyfert 2 galaxies, which we do not con- 
sider here. 

We found 7 Seyfert 1 galaxies with 0.01 < z < 0.07, 
listed in Table 3. The redshifts were remeasured in new 
spectra obtained with the ESO 1.52 m telescope, and cor- 
rected for heliocentric motion. Not all the fields from Ta- 
ble 1 were involved in this subsurvey, and the effective area 
was only 477 deg 2 (but see Sect. 4.2 for the computation 
of space densities). 

For these sources, proper photometry is at least as im- 
portant as for QSOs. The photometric systems of many 
earlier investigations of the luminosity function of Seyferts 
suffer from two severe drawbacks: (i) They were often 
based on heterogeneous collections of published (aperture) 
photometry, generally inconsistent with the magnitudes 
used for defining the survey flux limits, (ii) The contri- 
butions of host galaxies were either not corrected for at 
all, or some assumptions about universal intrinsic proper- 



ties of hosts and nuclei had to be made, such as colours 
(Sandage 1973) or host luminosities (Cheng et al. 1985). 

On the assumption that Seyfert 1 galaxies are the low- 
luminosity equivalent of quasars, the quantity of interest 
is the nuclear brightness, with the host galaxy contribu- 
tion removed. Although our small-aperture magnitudes 
were already much less affected by galaxy contamination 
than large- aperture or total brightness measurements, a 
major contribution could still be expected especially for 
the lowest nuclear luminosities. We therefore estimated 
these corrections individually, based on B and V images 
of the seven sample objects taken with the ESO 90 cm 
telescope. We first constructed empirical two-dimensional 
point-spread functions from nearby stars. These were then 
subtracted from the Seyfert images, scaled by a factor 
chosen such that the residual galaxy surface brightness 
distribution did not decrease inwards. The resulting host 
galaxy brightness was finally integrated over the aperture 
and subtracted from the total aperture flux. As can be 
seen in Table 3, the corrections were not negligible, and - 
expectedly largest for the lowest redshift sources. 

3. The surface density of bright QSOs 

An important diagnostic tool for QSO samples is the em- 
pirical relation between source fluxes and the number of 
sources per unit solid angle. Once the field properties and 
the sample photometry are complete, the cumulative sur- 
face density of QSOs brighter than magnitude B, N(< B), 
can be easily computed from summing over l/Sl c s f° r 
all relevant sample objects. From our QSO sample with 
z < 2.2 (47 objects) we have derived the A^(< B) re- 
lation shown as thick line in Fig. 2. Note that the ab- 
scissa values in this diagram are extinction-corrected mag- 
nitudes. The bright end of the curve is dominated by the 
extremely luminous source HE 1029—1401, but the main 
part can be well approximated by a simple power law of 
slope = 0.67 ± 0.04 (log N = /3B + const). If only the 
33 sources with z > 0.2 are considered, excluding thereby 
all QSO/Seyfert borderline cases (thin line in Fig. 2), the 
normalisation drops slightly but the slope is essentially 
unchanged with — 0.69 ± 0.05. For those readers who 
prefer differential surface densities, we have computed the 
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Fig. 2. Cumulative surface densities of bright QSOs with 
2 < 2.2. Thick line: this work, for z > 0.07, thin line: the same 
for z > 0.2. Open squares show the relation from Schmidt & 
Green (1983), modified for z > 0.07 (see text). Triangles give 
the LBQS (Hewett et al. 1995) relation, valid for z > 0.2. 

Table 4. Differential surface density A(B) of 33 HES QSOs 
with redshifts 0.2 < z < 2.2, expressed as number of QSOs per 
deg 2 per unit magnitude interval, n gives the actually observed 
number per half magnitude bin (B ± 0.25). 
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figures in half magnitude bins between B = 14.75 and 
B = 17.25 (Table 4). The counts were multiplied by a 
factor of 2 to formally express number counts per unit 
magnitude intervals. 

Comparison relations from other surveys are numerous 
for fainter magnitudes, but rare for B 16.5. The Palo- 
mar Bright Quasar Survey (BQS; Schmidt & Green 1983, 
hereafter quoted as SG83) provides the only available op- 
tically selected QSO sample with a well-defined flux limit 
subtending over a similar region in the Hubble diagram. 
SG83 listed differential surface densities computed from 
the full set of 114 QSOs with z < 2.2 found over an area 
of 10 714 deg 2 . To make their numbers comparable to ours, 
we recomputed the BQS surface densities by removing all 
z < 0.07 entries from their sample, and converted the re- 
sult into a cumulative relation. This relation is shown in 
Fig. 2 by the open squares. A large discrepancy at all mag- 
nitudes is apparent. At B = 16, the surface density found 
in the HES is high by a factor of 3.6 compared to the 
BQS, certainly more than can be permitted by statistical 
fluctuations. 

A similar result was already reported by Goldschmidt 
et al. (1992) from an analysis of the Edinburgh sur- 



vey, who found 5 QSOs above the BQS flux limit in 
an area of 330 deg 2 (non-overlapping with the present 
HES area) while the BQS contained only one. Their es- 
timate of N{B < 16.5) = 0.024 deg- 2 for 0.3 < z < 2.2 
is completely consistent with our corresponding value of 
0.021±0.05. Furthermore, Goldschmidt et al. found a zero- 
point offset of 0.28 mag between the BQS photographic 
magnitudes and theirs, in the sense that the Palomar mea- 
surements were systematically too bright, thus increasing 
rather than relaxing the discrepancy. 

It should be noted that neither SG83 nor Goldschmidt 
et al. corrected their magnitudes for Galactic extinction, 
whereas we have done so. However, our survey area is pre- 
dominantly located in regions of rather high As, while 
SG83 estimated that the extinction for most BQS quasars 
were 'not much more than 0.1 mag'; this even might be, 
at least partly, compensated by the above mentioned zero- 
point offset. Because of these uncertainties we have not at- 
tempted to correct the BQS results for extinction. One can 
estimate globally that an average j4b(BQS) of 0.61 mag 
would be necessary to bring BQS and HES surface densi- 
ties to a match - a totally unrealistic value. 

At fainter magnitudes the agreement of HES num- 
ber counts with those of other surveys is excellent. As 
an example, we have plotted the values from the LBQS 
(Hewett et al. 1995) into Fig. 2, with the abscissa shifted 
by +0.1 mag to transform their Bj photometric system 
approximately into B magnitudes. While the counts in 
the brightest bins are again lower than ours, this can be 
understood as a result of poor statistics in the LBQS for 
B ^ 16.5, possibly arising from saturation effects. Around 
B ~ 17, HES and LBQS surface densities join smoothly, 
without any detectable significant offset. Similarly good 
agreement is reached with other surveys, cf. the recent 
compilation by Cristiani et al. 1995. We are thus confident 
that the photometric scale used in the present investiga- 
tion is essentially unbiased and adequate for luminosity 
function work. 

4. Luminosity functions 

4-1. The luminosity function of low-redshift QSOs 

For all objects in Tabic 2, we first estimated absolute mag- 
nitudes for the rest-frame B band by computing their lu- 
minosity distances in an expanding Fricdmann universe 
with the formula of Terrell (1977), and applying K cor- 
rections taken from Cristiani & Vio (1990). As low-redshift 
subsample we selected the 20 QSOs with z < 0.3. Space 
densities were derived using the generalised 1/V max esti- 
mator (Felten 1976; Avni & Bahcall 1980), thus incorpo- 
rating the full information about field-dependent flux lim- 
its, including fields with zero detections (the 'coherent' 
analysis in the terminology of Avni & Bahcall). 

The resulting cumulative luminosity function is shown 
in Fig. 3. Locations of individual objects are indicated by 
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Fig. 3. Cumulative luminosity function of QSOs with 
0.07 < z < 0.3 (thick line). For comparison: Cumulative lu- 
minosity function of BQS quasars with 0.07 < z < 0.3, con- 
structed as described in the text (dotted line and small sym- 
bols). 
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Fig. 4. Cumulative luminosity function of the additional 
Seyfert 1 sample (z < 0.07). Thick line: LF for nuclear mag- 
nitudes corrected for the host galaxy contributions; thin line: 
uncorrected LF. For comparison: Seyfert 1 nuclear luminosity 
function of Cheng et al. (1985; dotted line and small symbols) 



the error bars. We have chosen the cumulative represen- 
tation ip(< Mb) as it allows one to avoid binning (par- 
ticularly problematic for smaller samples) and to make 
the contribution of each object in the sample apparent. 
The errors were estimated from Poisson statistics in the 
sample, combined with the uncertainties in volume deter- 
mination resulting from the photometric errors. 

Restricting the HES luminosity function to z < 0.2 (14 
objects) results in slightly larger error bars, but the space 
densities change only insignificantly, by less than 0.05 dex, 
and, ovcrplotted in Fig. 3, could not even be distinguished 
from the z < 0.3 relation. We thus feel justified to neglect 
differential evolution within our rcdshift shell; we shall 
further investigate and qualify this assumption below. 



Up to now, the BQS provides the only comparison sam- 
ple in the region Mb < —23 and z < 0.3, containing 45 
objects that fulfill also the condition z > 0.07. Surpris- 
ingly, there is yet - to our knowledge - no direct deter- 
mination of the local QLF based on BQS objects alone. 
We have therefore computed a QLF from this sample with 
the same methods used for our own data, taking magni- 
tudes, survey limits, and areas from SG83. We show the 
results of these computations as a cumulative distribution 
in Fig. 3. (Since all other literature data are available only 
in binned form, we have also binned the BQS LF in Mb, 
to achieve a more uniform presentation.) As mentioned 
above, Galactic extinction is neglected for these data; at 
any rate, the effect should be very small. 

The space densities found from the HES data are much 
higher than those derived from the BQS, with the discrep- 
ancy increasing with luminosity. For the most luminous 
sources, we find that there are almost an order of magni- 
tude more low-redshift quasars per unit volume than dis- 
covered by the BQS. Possible reasons for this discrepancy 
are discussed in Sect. 5. 

4-2. The space density of Seyfert 1 nuclei 

To extend the local luminosity function towards fainter 
levels of nuclear activity, we have used the additional 
Seyfert 1 sample with z < 0.07. The absolute magnitudes 
were computed in the approximation of Euclidean static 
geometry and neglecting K corrections. To obtain reli- 
able and unbiased estimates for the accessible volumes, 
the maximum redshift z maXji at which object i would still 
have been included in the survey had to be determined for 
each object. Because of the superposition of a host galaxy 
with constant intrinsic scale length and a point-like nu- 
cleus, the small-aperture magnitudes depended on redshift 
in a complicated way, and z maXj i was therefore not only a 
function of B and Bn m but also of intrinsic properties. We 
used the available CCD images to simulate the mixing of 
AGN / host contributions into the aperture as a function 
of rcdshift z, resulting in an array of z max ,fc values, one 
for each field k. For simplicity, Table 3 supplies 'effective' 
values z max that permit a direct conversion into space den- 
sities. A V/Vmax test (Schmidt 1968), implemented in the 
generalisation proposed by Avni & Bahcall (1980), gives 
a mean value of V/V max = 0.43 ±0.11 (individual values 
are also listed in Table 3). Compared with the expectation 
value of 0.5 for a complete, non-evolving, homogeneously 
distributed sample, there is no evidence for major incom- 
pleteness. 

The resulting luminosity function is displayed in Fig. 4. 
It can be seen that a host galaxy correction is impor- 
tant particularly for the intrinsically faintest objects in 
the sample: Without correction, the luminosity function 
shows a steep upturn at Mb — —20; for these objects, 
even the small-aperture magnitude is dominated by the 
host rather than by the AGN. 
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Fig. 5. Combined local luminosity function of quasars and 
Seyfert 1 nuclei with z < 0.3, based on magnitudes corrected 
for the host galaxies. 



Table 5. Binned differential local luminosity function, com- 
puted from 27 QSOs and Seyferts with z < 0.3, in Mpc -3 per 
unit absolute magnitude interval centred on the value given by 
Mb- Column n gives the number of objects per bin. 
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Although the sample is not large, the derived space 
densities are fully consistent with previous estimates 
(Cheng et al. 1985; Huchra & Burg 1992). However, none 
of these previous investigations have payed as much at- 
tention to an appropriate host galaxy correction, neither 
for source luminosities nor for the survey volume, as we 
have. In particular, the analysis of the CfA sample by 
Huchra & Burg was based on uncorrected Zwicky mag- 
nitudes (essentially measuring the total galaxy bright- 
ness), making a straightforward comparison impossible. 
More similar to our approach was the work of Cheng et 
al. (1985) with Markarian galaxies, but their sample, al- 
though much larger than ours, had to be manipulated with 
substantial incompleteness corrections, in addition to the 
afore mentioned problems of a heterogeneous photometric 
data base. 

4-3. The combined local luminosity function 

The local universe is the only domain where it is possible 
to construct the entire QLF, including low-luminosity tail, 
directly from one single survey. The local luminosity func- 



Fig. 6. The bright end of the quasar luminosity function up 
to 2 = 0.5, in comparison with the local QLF. Shown are the 
cumulative relations in the redshift shells z < 0.2 (thin line), 
and 0.2 < z < 0.5 (thick line). 



tion of QSOs and Seyfert 1 nuclei, obtained by combining 
the two independently derived luminosity functions pre- 
sented above, is valid between absolute magnitudes —18 
and —26, thus spanning three decades in luminosity. 

Overlaying the cumulative relations of Figs. 3 and 4, 
the match in the common luminosity region is remark- 
ably good (cf. Fig. 5). To arrive at a consistent system of 
nuclear luminosities, the magnitudes in the QSO sample 
were corrected by subtracting a template host galaxy of 
Mb = —21. Since the magnitudes had already been mea- 
sured through a small aperture, this host subtraction was 
a minor correction for all objects in the QSO sample, and 
the results are not sensitive to the assumed host luminos- 
ity. 

The combined local QLF can be fitted well by a 
simple linear relation between Mb and log^(M), corre- 
sponding to a single power-law for the differential QLF 
{4>{L)dL oc L a ). A slope parameter a = —2.17 ±0.06 pro- 
vides a statistically acceptable fit over the full range of 8 
magnitudes. Fitting only the QSO sample (z > 0.07) does 
not alter the slope significantly (a = — 2.10±0.08); fitting 
only the z < 0.07 Seyferts yields a = —2.38 ± 0.23, again 
consistent with a constant slope of a ~ —2.2. The slope 
of the BQS LF is much steeper: Apart from flattening at 
both high and low luminosities, the data demand a Ss —3 
over the most of the range. 

From the combined dataset, we have also produced a 
binned differential luminosity function 4>(M), tabulated 
in steps of one in absolute magnitude (Table 5). As was 
to be expected, this representation is quite sensitive to 
the actual choice of binning intervals, most obvious in the 
succession of a 'high' and a 'low' bin at M B = —22.5 and 
—21.5, respectively. There is, however, full consistency of 
the results obtained; in particular, the QLF slope param- 
eter a ~ 2.2 is reproduced within the expected uncertain- 
ties. 
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4-4- The quasar luminosity function at higher redshifts 

The sample of Table 2 contains eight QSOs with redshifts 
0.3 < z < 0.5 that we used to further constrain the poten- 
tial error introduced by neglecting differential evolution. 
Binning the full z < 0.5 sample into two redshift regimes 
0.07 < z < 0.2 and z > 0.2, respectively, there are 14 ob- 
jects in each shell. The luminosity functions are compared 
in Fig. 6: For -25 <> M B & -24, the local and the ad- 
jacent 0.2 < z < 0.5 shell have essentially identical space 
densities, while for Mb < —25, the QLF in z > 0.2 seems 
to be merely a continuation of the local QLF to higher 
luminosities, with the same slope and normalisation. Di- 
viding the sample at z — 0.3 instead of 0.2 does not affect 
the results. 

Outside the local universe, a single survey can yield 
no more than a small segment of the QLF, and only the 
combination of many different surveys can provide simi- 
lar coverage of source luminosities as in the local QLF. 
The HES has the potential to provide new samples cov- 
ering the brightest parts of the known QSO population, 
and it will be a promising task for the future to combine 
our results with those of fainter surveys. Such an analysis 
is, however, not the scope of the present paper. We only 
remark that even at redshifts higher than z = 0.5, the 
QLF pieces sampled by the HES do not show significant 
deviations from an extrapolation of the local relation. A 
more quantitative discussion of this issue is deferred to a 
later paper. 

5. Discussion 

5.1. Incompleteness of the BQS 

Our results imply that the Palomar Bright Quasar Sur- 
vey is not only incomplete, but also heavily biased. The 
general incompleteness may be estimated from the sur- 
face densities to be around a factor of 2-3. This effect 
has been noted by others before (Wampler & Ponz 1985; 
Goldschmidt et al. 1992), and we have at present no easy 
explanation for it. Some mild adjustment of the BQS data 
to higher values of space densities and luminosities may be 
required because of neglected Galactic extinction, but this 
cannot account for the observed discrepany. One possibil- 
ity is that the rather large photometric errors in the BQS, 
reflected in the (U — B) colours, scattered many QSOs out 
of the UV excess domain. 

There seems to be an additional deficit of luminous 
low- redshift QSOs, by another factor of > 2 for Mb < 
—24. This discrepancy is apparent already from simply 
counting the number of z < 0.3 QSOs more luminous than 
Mb = —24.5: There are 5 objects in the present HES sam- 
ple, while the entire BQS contains only 15 such sources, 
in a more than 15 times larger area. As this is the first 
detection of such an effect, one has to carefully consider 
systematic errors. Note, however, that the HES is the first 
optical QSO survey since completion of the BQS that is 



sensitive to a similar range of redshifts and luminosities, 
and it is therefore not altogether inconceivable that sub- 
stantial selection biases have remained unnoticed over a 
long time. 

The good agreement between the HES z > 0.2 or 
z > 0.3 surface densities with other major surveys, e.g., 
the LBQS (Hcwett et al. 1995), the Edinburgh survey 
(Goldschmidt et al. 1992), or the HBQS (Cristiani et al. 
1995), makes it implausible that the HES photometric ze- 
ropoint could be seriously in error. 

We also tested our method of zero-point corrected 
CCD photometry All relevant quantities were recomputed 
based on the standard approach of using only photo- 
graphic magnitudes measured over large apertures, with- 
out any systematic change of results, except that the lu- 
minosities of low-rcdshift QSOs with bright host galaxies 
were often overestimated in the photographic data, cre- 
ating a bias towards even higher QLF values. The CCD 
approach thus seems justified and appropriate. 

The reconstruction of the 0.07 < z < 0.3 luminosity 
function for the BQS was based on the same procedures as 
applied to our own data. We neglected differential evolu- 
tion within the z < 0.3 shell, but the observed discrepancy 
is independent of any assumed evolution law, and there- 
fore most likely caused by real incompleteness in the BQS 
sample. 

The exclusion of objects of non-stellar appearance in 
the BQS is an intuitively appealing candidate mechanism 
for this additional incompleteness. It is obvious at first 
glance (and has in fact never been disputed) that the BQS 
sample is systematically lacking low-luminosity Seyferts; it 
has usually been assumed that this selection effect ceased 
to be valid for Mb — 23. Our results suggest that this 
is not the case. However, it is then somewhat surprising 
that the incompleteness should be smallest, almost dis- 
appearing, in the luminosity range around Mb — —23, as 
seems to follow from Fig. 3. In fact, this would then be the 
only region of the entire Hubble diagram where the BQS 
counts were not significantly below those of the HES. 

Increased incompleteness for the most luminous QSOs 
seems a possible, but rather unlikely explanation. If, in- 
stead, the BQS magnitudes particularly for the slightly re- 
solved Mb — —23 objects were systematically too bright, 
overcompleteness would occur, scattering objects into the 
sample from below the flux limit. We know from our own 
experience with photometry of resolved objects on digi- 
tised photographic plates that such effects exist, especially 
if a large digitisation aperture is used (this is why we de- 
cided to incorporate the CCD measurements). The BQS 
was based on PDS scans with an aperture of 4'.'4 (Green 
& Morrill 1978), and magnitudes were basically isopho- 
tal in photographic density. We suspect that these techni- 
cal details provide a clue to the above contradiction, and 
that incompleteness due to the rejection of extended QSO 
hosts could be partly counteracted by photometric bias 
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and overcompleteness for only marginally resolved objects. 
Further analysis of this point is certainly required. 

If our hypothesis is correct, two important conjectures 
can be made: (1) The slope of the QLF for M B < -23 is 
much flatter than previous investigators estimated, weak- 
ening the evidence for a 'break' in the Seyfert/QSO tran- 
sition regime (e.g., Marshall 1987) - the local QLF pre- 
sented here showing essentially no break at all. Note that 
also in the Seyfcrt LF of Cheng et al. (1985), the pres- 
ence of such a break is not really demanded by the data. 
(2) The BQS sample, up to now the prime source for lu- 
minous radio-quiet low-redshift QSOs, might suffer from 
severe selection biases with respect to the morphological 
properties of their galaxy hosts. 

5.2. Constraints on quasar evolution 

The local luminosity function derived in this paper is in- 
dependent of samples collected at higher rcdshifts, and 
therefore well suited to test and constrain possible evolu- 
tion laws. Existing estimates of the QLF at rcdshift z ~ 2 
(e.g., Boyle et al. 1991) show a marked break, or change 
of slope, at intermediate luminosities, that shifts towards 
lower values with decreasing z. We have argued above that 
the local equivalent of this break is probably an artefact. 
The flat slope of the z « QLF derived in this paper, and 
the absence of clear features, demand, on the assumption 
that the high z QLF determinations are correct, that both 
shape and slope of the QLF change strongly with redshift. 
This clearly is in contradiction to the 'standard picture' of 
pure luminosity evolution (PLE), where a double power- 
law luminosity function with invariant logarithmic shape 
is merely shifted along the luminosity axis. 

However, it should be noted that virtually all evolu- 
tion modelling since 1983 had to rely on the BQS data of 
SG83, to cover an otherwise inaccessible part of the Hub- 
ble diagram. Hence, in these regions, the models merely 
reproduce the incompleteness and possible biases of the 
BQS contribution. The problem that QSO surveys were 
notoriously incomplete at low redshifts was of course re- 
alised by most researchers, and it has become customary 
to apply a lower redshift cutoff of z min ~ 0.3 to all samples. 
The construction of a z = QLF from such models has 
therefore always the character of an extrapolation outside 
the validated range, and should be treated with appropri- 
ate caution. Even so, a quantitative comparison between 
evolutionary models and actual measurements of the lo- 
cal Seyfert 1/QLF is a powerful test which deserves higher 
attention than it has been paid to in the past, compared 
to the careful analyses of higher redshift data. 

In the construction of the z < 0.3 local QLF we have 
not corrected individual luminosities for statistical evolu- 
tion. While it is true that such a neglect could artificially 
flatten the luminosity function, evolution must be fairly 
strong for the effect to be relevant. We show now that the 
no-evolution assumption was, in fact, the best approxi- 



mation one could possibly make, over the range of lumi- 
nosities covered by the HES. In the case of evolution, the 
LFs of adjacent, non-overlapping redshift shells should be 
different. For the simple case of PLE with the rate sug- 
gested, e.g., by Boyle et al. (1991), one would expect an 
abscissa offset of AM = —0.8 mag between z < 0.2 and 
0.2 < z < 0.5. The necessity to invoke such an offset is not 
apparent from Fig. 6. While it might be consistent with 
our data to allow for a somewhat steeper intrinsic slope of 
the QLF, combined with mild evolution, there is certainly 
no evidence for evolution in the data. 

These conclusions are supported by results reported 
recently from other surveys. Hewett et al. (1993), with a 
preliminary analysis of the LBQS, showed that the QLF 
flattens considerably at z < 1; a similar conclusion was 
reached by Miller et al. (1993) using the Edinburgh sur- 
vey. Although both surveys have the usual lower redshift 
limit excluding the local population (z > 0.3 for the Edin- 
burgh survey, z > 0.2 for the LBQS), a prediction for the 
slope of the local QLF can be made by extrapolating the 
trend seen by Miller et al.; the result is consistent with 
our measured value of a = —2.2. A flatter QLF implies a 
higher fraction of luminous quasars among the total pop- 
ulation. Our results indicate therefore, in agreement with 
Hewett et al. and Miller et al., that the evolution of the 
most luminous quasars must proceed considerably slower 
than suggested by the notion of pure luminosity evolution. 

6. Conclusions and outlook 

We have analysed a new, flux-limited and well-defined 
sample of QSOs & Seyfert 1 galaxies. Host-galaxy depen- 
dent selection and photometric biases, if not absent, are 
greatly reduced in comparison to other optically selected 
samples. We find that the derived space density of lu- 
minous QSOs in the local universe is much higher than 
previous surveys indicated, and that consequently these 
objects show a much slower cosmological evolution. 

With the present sample of 55 QSOs and Seyferts in 33 
fields, only a small fraction of the HES area is covered. We 
have already acquired, and partly processed, much addi- 
tional plate material enlarging the survey by a significant 
factor. For example, we now have several fields in common 
with the BQS, and it shall soon be possible to pursue the 
question of incompleteness and biases in the BQS by a 
direct comparison of objects and photometry. 

To understand and model the process of quasar evo- 
lution requires strong constraints from observations. It 
seems now that we are farther away from a coherent em- 
pirical picture of the QSO population and its evolution 
than thought a few years ago. Quasars of different lumi- 
nosities evolve differently, and quasars do not necessar- 
ily evolve with the same rate at different cosmological 
epochs. There are presumably more parameters needed 
to describe quasar evolution properly than in the simple 
model of 'pure luminosity evolution'. 
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The ultimate aim of QSO survey work is the deduction This article was processed by the author using Springer- Verlag 
of intrinsic physical properties from statistical samples. MgX A&A style file L-AA version 3. 
Unfortunately the physical processes in active galaxy nu- 
clei are poorly understood, and the relations between sta- 
tistical and physical properties therefore far from unique. 
The best study cases for nuclear activity in galaxies are 
nearby Seyferts; as shown above, the luminosity function 
of Scyfert 1 nuclei can be smoothly continued into the clas- 
sical quasar regime, without a significant change of slope 
or break. While this by no means proves that Seyfert 1 
nuclei are simply scaled-down versions of quasars, it sug- 
gest a continuity of properties rather than distinct classes. 
With the completion of the HES, a large and well-defined 
sample of highly luminous low-redshift QSOs will become 
available, to study the relationship between QSOs, their 
environments, and their evolution in detail. 
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